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Abstract. Spline wavelet tight frames of [20] have been used widely 
in frame based image analysis and restorations (see, e.g. survey ar- 
ticles [15ll22j ). However, except for the tight frame property and the 
approximation order of the truncated series (see |13ll20p . there are few 
other properties of this family of spline wavelet tight frames to be known. 
This paper is to present a few new properties of this family that will 
provide further understanding of it and, hopefully, give some indications 
why it is efficient in image analysis and restorations. In particular, we 
present a recurrence formula of computing generators of higher order 
spline wavelet tight frames from the lower order ones. We also represent 
each generator of spline wavelet tight frames as certain order of deriva- 
tive of some univariate box spline. With this, we further show that each 
generator of sufficiently high order spline wavelet tight frames is close to 
a right order of derivative of a properly scaled Gaussian function. This 
leads to the result that the wavelet system generated by a finitely many 
consecutive derivatives of a properly scaled Gaussian function forms a 
frame whose frame bounds can be almost tight. 



1. Introduction 

This paper is to investigate the family of the sphne wavelet tight frames 
derived from [2Uj . We start with basic notions. For given ^ := {^i, . . . , -0^} C 
L2(M), the wavelet system generated by "if is defined as 

X{^) := {iPe^^^k := 2"/V£(2" ■ -k) : 1 < i < r; n, k e Z}. 
The system X{^) C L2(M) is called a tight frame if 

/= ^ if, 9)9 

holds for all / G L2(M). If C L2(K) is a tight frame system of L2(M) 

generated by a Multiresolution analysis (MRA), then its generators ^ are 
called framelet. 
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The MRA starts from a refinable function (p. A compactly supported 
function (p is refinable if it satisfies a refinement equation 

(1) (^(x) = 2j^a,v9(2x-j), 

for some sequence a S ^2(^)- Refinable equation ([T|) can be written via its 
Fourier transform as 

^(^^) = a(^) • a.e. w G M. 

We call the sequence a the refinement mask of (/? and a{(jj) the refinement 
symbol of ip. Here, we use / to denote the Fourier transform of / G -Li(M), 
which is defined as 

/oo 
f{x) e'Kp{—iojx)dx. 

For a refinable function ip G L2{^), let Vq be the closed shift invariant 
space generated by {ip{- — k) : k £ Z} and Vj := {f{2^-) : f G Vq}, j G Z. 
It is known that when p is compactly supported, then {Vj}j^z forms a 
multiresolution analysis. Recall that a multiresolution analysis is a family 
of closed subspaces {Vj}j^z of -^2(1^) that satisfies: (i) Vj C V^+i) (ii) Uj 
is dense in L2(M) and (iii) flj Vj = {0} (see [2] and [H]). 

A special family of refinable functions is B-splines. Let yj^™') be the cen- 
tered B-spline of order m, which is defined in Fourier domain by 

(2) ^MM = e-^sinc(|r, 
where 

I sin(x)/x, for X 7^ I 0, m is even 

(3) smc(x) := < ; and Jm '■= \ 

I 1, for X = I 1, m is odd 

Then is refinable with refinement symbol 

aM(a;)=e-^cos-(|). 

The tight framelets can be constructed by the unitary extension principle 
(UEP) of [20] from a given multiresolution analysis. For a given B-spline 

ipi^) of order m, it was shown in |20j that the m functions, ^f^™) = {ip^"^^ : 
i = 1 , . . . , m} , defined in Fourier domain by 

ir(^) - ,v-j-,M °"'""'t):;r'"^^' . 




form a tight wavelet frame in L2(M), i.e. ^('"^ is a framelet set. We call 
^i^"^^ as the B-spline framelet of order m. The B-spline framelet is either 
symmetric or anti-symmetric and has small supports for a given smoothness 
order. Similar with B-splines, each B-spline framelet has an analytic form. 
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Since the publication of the unitary extension principle (UEP) of [20] 
in 1997, there are many theoretic developments and applications of MRA 
based wavelet frames. In particular, the B-spline framelets ^f^™-) derived 
from the UEP in [20] are widely used in various applications, which include 
image inpainting in [5]; image denoising in [8]; high and super resolution 
image reconstruction in [10]; deblurring and blind debluring in [HHS]; and 
image segmentation in [H] . The interested reader should consult the survey 
articles [l5l[22] for details. 

The paper is organized as follows: Section 2, which contains two subsec- 
tions, is for some basic properties of B-spline framelets. In particular, in 
sub-section 2.1, we present recurrence formulas for B-spline framelets 
in which the well-known recurrence formula of B-splines can be viewed as a 
special form of recurrence formulas of B-spline framelets. This gives a fast 
algorithm for computing them. We further show that the B-spline framelets 
can be derived from the ith. derivative of some univariate box spline in sub- 
section 2.2, This was implicity used in [4], to approximate the some deriva- 
tives of a function. In Section 3, we investigate the asymptotic property of 
B-spline framlets ip£^\ i = 1, . . . ,m. We firstly prove that the univariate 
box splines defined in Section 2 uniformly converge to a Gaussian function 
under a mild condition, and we further show that 

max max|^l™^(x) - G^J^\x)\ < ,\ — , 
i<e<m xm *^ * \ /I ~ ^3/2 

where G^"^^ is the £th derivative of some scaled Gaussian function G{x) (see 
Section 3.2 for the detailed definition). 

This leads to discover that wavelet system generated by a finite number 
of consecutive directives of scaled Gaussian function form a frame whose 
bounds are almost tight, and that is done in Section 4. 



2. Properties of B-spline framelets 

In this section, we give a recurrence formula of the B-spline framelets 
which computes higher order framelets from lower order ones. We also 
show that one can represent the derivatives of higher order framelets by 
lower order framelets. Furthermore, we derive another set of formulas that 
represents each framelet as a derivative of a univariate box spline which is 
already implicitly used in [4], where a theory is developed to connect the 
PDE based and spline wavelet based image restorations. 

2.1. Recurrence formulas of B-spline framelets. While the recurrence 
formulas for B-splines and their derivatives are well-known (see [1]), the 
corresponding formulas for B-spline framelets are not available yet. This 
section is to establish such formulas. Let Bm ■= 93^™"^ (• -|- jm/2), where 
(^(™-) ig given in ([2|) and jm is defined ([3]). Recall the following well-known 
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recurrence formula of B-splines: 

,, 2x + m + l f l\ m + l- 2x f 1 
(5) Bm+i[x) = Bm. a; + - H Em [x- - 



2m V 2y 2m \ 2 

Based on ([5D, one can compute B-splines fast and easily which makes B- 
splines useful. The derivative of B-splines can be computed by the lower 
order splines as given below: 

(6) ^S„,+i(x) = B.m + ^) - ^™ - ^ 

The aim of this section is to give the corresponding formulas for the B-spline 
framelets V'^™^ ,£ = 1, . . . ,m. To state the formulas conveniently, we present 
the formulas for the function := + Note that the Fourier 

transform of ipi"^^ is 



^^^^cos™-^(|)sin-+^(f) 
(7) (w) = 2 ■ " ' 




m 



We note that the formulas presented in this subsection are used to calculate 
the function value and the derivative of V"!™^ • When m is even, V'^"^^ = ^^P^ ■ 
When m is odd, one can obtain those of the function ^/'^'"^ by the half- 
translation of V'^™^. Hence, the formulas given in this subsection also work 
for V'^™^ with a proper shift. 

Next, we present the recurrence relations of framelets 

Theorem 1. Let m N be given and 1 < £ < m and let the framelet ip^"^^ 
derived from B-spline of order m be given via its Fourier transform as 



Then, we have the following recurrence formula between ip^"^'^^^ and ipf^ : 
for 1 < i < m: 



771 + 1 

m + l-l 



2x + m+l ~Un) ( 1 \ m + 1 - 2x / 1 \ ^ lirn), . 

i)) s + - H ip) 2; - - H W) (x) 



2m * V 2y 2m *^ V 2J m 

the recurrence formula between 'ip'^J^j^^^ and tpin^^ is: 
(9) 



Proof. We firstly prove ^ which is done in Fourier domain. Note that 

d ■ 
duj 
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which imphes that the Fourier transform of function ge{x) := xt/j^^^\x) is 
(10) 

Note that 



2?Ti ' 2' m"^ 2' 2"'^ ' 2' 

and 

Zm 2 2 2 m 2 

A simple manipulation shows that the Fourier transform of the right hand 
side of dSl) becomes 



771 + 1 ,1 , , 1 , .iuj. . iuj..^("^). . 

-,{— exp(— j^Hw) + - • (exp(— ) + exp{-—))i;^ [u) 



m+l-tm 2 ' ' ' 2 ' '2' ' 2 

— exp(-— + —i)^ (w)) 
m 2 m 



,/m + l\cos"^+i"^(^)sin"+i+^(^) -{rn+i) 



This proves ([8]). Similarly, the Fourier transform of the right side of ([9]) is 

— ( exp(y)(5™(a;) + y^/;^ (cj)) + exp (-y - y^/;^ (uj)) - 2gm{uj) 

which proves ([9]). □ 

Furthermore, combining ([8]) and ([9]), we have a recurrence algorithm for 
efficiently computing = 1, . . . ,m. When £ < m, we can use 1^ to 

compute ijjf^^ by ip^^ we can use ([9]) to compute ipf^ by ipf_-^\ Hence, 
we finally can reduce the computation of ip£^^ to that of "^^^ . Note that the 
function ip^^ is Haar wavelet and 



1, if [-1/2,0) , 
-1, if x G [0,1/2] , 
0, if |x| > 1/2. 



We next show the method for computing ijj^^ by a table. In the following 
table, for the notation — )•, we use the formula ([8|), while for the notation \, 
we use (l9l): 
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Bi 


B2 


B3 


Ba 










\ 


- 







Using the method, we compute B^ and corresponding framelets (See Fig- 
ure 1). 




50 100 150 ' 50 100 150 ' ED 100 150 




Figure 1. The and corresponding framelets. 

Next, we give the recurrence formula for the derivatives of '>p£^^- 

Theorem 2. Let m £ N be given and 1 < £ < m, and let the framelet ip^"^^ 
derived from B-spline of order m be defined by its Fourier transform as 
When 1 < i < m — 1, we have 



d ~{m), ^ _ / rn f j{m-l) ( 1\ 7(m-i) / 1 



When £ = m, we have 



Proof. We prove (lllh here while (I12p can be proven similarly. A simple 
calculation shows that the Fourier transform of the right side of pip is 



(13) 




cos'"-^-^(|)sin'^-i+^(|) 



(e 2 — e 2 



.4, 



(I) Sin 



(f) 



m— 1 
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Note that the Fourier transform of is 

Combining (jl3p and (jl4p . we conclude (jlip . □ 

Remark 1. Note that -(/i^""'' = Bm- If we take £ = in ([8]), the recurrence 
relation ([8]) is reduced to ([5]), which is the recurrence formula for B-splines. 
Similarly, if we take £ = in (fTT]l . then (fTTI) is reduced to the derivative 
formula of B-splines (j6]) . 

2.2. Representing ^Z'^™^ as the ^th derivative of a univariate box 

spline. We first recall the definition of box splines. The hox spline B{-\'S) 
associated with a matrix H G M'^^" is the distribution given by the rule 
(see m) 

(15) / B{x\'S)ip{x)dx = I ip{Eu)du, for all if G &{R'), 

where is the test function space. The box spline can be consider as 

a volume function of the section of unit cubes (see [3ll241[25]). If we take 
S = (1, 1, . . . , 1) G M-^^*", then the box spline B{-\E) is reduced to a B-spline 
of order m. In the following theorem, we show that the B-spline framelet 
can be considered as the higher order derivative of a box spline up to a 
constant. 

Theorem 3. Let m G N 6e given and 1 < i < m. Suppose that the framelet 
'(/'^'"^ is defined by its Fourier transform in Qj. Set 




Then 

(16) ^^^(-)=/(T)-i-|z^(--fi^-^)' 

where jm is defined in In particular, ipln^^ is the m-order derivative of 
-S2m(2 • — jm)/4"*, where B2m is the B-spline of order 2m. 

Proof This again is proven in Fourier domain. It follows from the definition 
of box splines (|15p that the Fourier transform of the box spline B^-lEm/) is: 



B{u}\Em,e) = 
Then the Fourier transform of 



sm^\™-^ /sinf ^ 2f 
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can be computed as: 




21 




which proves (jl6p . According to the definition of box sphnes, we have 



B2m{2x - jm) = B{x - 



And hence, V'm^^ is the m-order derivative of i?2»n(2 • — Jm)/4'". □ 

Theorem [3] states that each framelet of vanishing moment order of i is 
the ith derivative of a univariate box sphne whose support is the same as 
the framelet and whose Fourier transform dose not vanish at origin. This 
fact is used in |3j to discretize differential operators by using framelets. We 
illustrate here for the case where m is even how the framelet coefficients 
can be viewed as the samples of a difference of a given function which can 
be used to approximate the derivative of this function when it is smooth. 
Recall that, if / G L2(M), then 

m 

i=l k,n€Z 

We next show that the coefficients (/, ipe^k) is the ^-order difference of a 
discretization of / up to a constant. We set 

Tnfix) :=2t^/(^),i?(.|H:„,,) 

where 

-=11 1 i h 

^m,£ ■ l ^^ • J ' y 2 ' ' ' ' ' 2 

m-i ' ^ ' 

Then T„/ can be considered as a discretization of /. We define the difference 
operator by 

Air„/(x) := T„/ U+l)- Tnf (x--^ 



ON B-SPLINE FRAMELETS DERIVED FROM UEP 



9 



Now we can give an explanation of the coefficient of (/, differ- 
ence operator Ai. In fact, we liave 

4 

4£ ^ \ V 2" / 

4^ yhi/2,i/2)-+. V 2" ; 

Here, the first identity follows from the result of Theorem [3] and the inte- 
gration by parts. The second relation is obtained by the definition of box 
splines. 

Remark 2. Theorem [3] shows that one can obtain the B-spline framelet 
by calculating the derivative of box splines, which provides a new path 
to construct spline framelets. We hope to construct multivariate spline 
framelets by calculating the derivative of some relevant box splines in future 
work. 




3. The asymptotic property of B-spline framelets 

3.1. The asymptotic convergence of univariate box splines. Up to 

the normalization, the B-spline tends to Gaussian function pointwise and 
in all Lp norms with 2 < p < +oo as the order tends to infinity (see |23j). 
Motivated by the results, in this subsection, we investigate the asymptotic 
convergence of univariate box splines, which is helpful to understand the 
convergence of , with V'^™^ being the £-order derivative of a box spline 
up to a constant. 

To state the results conveniently, throughout the rest of this paper, we 
shall use the notation X <a,6,... Y to refer to the inequality X < C -Y , where 
the constant C may depend on a, b, ... , but no other variable. In the next 
theorem, we show that the normalized box splines converge uniformly to a 
Gaussian function. 

Theorem 4. For each /c G N, Let 

be a given a set of points with a^-^^ > 0, , j = 1, . . . , A;. Let B{-\Ek) be the box 
spline associate with H^. Assume that 

(17) \\Ek\\l = a"^ + €k, 
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with a is a fixed constant and limfc_j>oo ek = 0, o,nd assume that 

(k) 

(18) Ci < < C2 

mini<j<fc 

where c\ and C2 are fixed positive constants independent of k. Then, 

(19) 

6 f 6x^\ ^, ,„ , I . {Inkf 



In order to prove TheoremSl we need the following technical lemma about 
the Fourier transform of the box spline i?(-|Hfc): 

Lemma 5. Under the conditions of Theorem^ 

, , , / / (lnA;)2 

(20) max|/fc(a;) - exp ( — I | <ci,c2 — ^ ^ Kfcl • |ln|efc||, 



where 



^ ( 
fk{oj) ■■= B{u)\Ek) = JJsinc 



(k) ^ 
a) 'iv 



Proof. Without loss of generality, we suppose that for each fixed k 

0<af)<af)<...<af. 
Then ^ and ([H]) imply that 

We firstly consider the case > n /a^^\ Note that sinc(-) is a monotone 
decreasing function in [0, vr] and 

|sinc(cij)| < — , for |a;| > vr. 

Then, we have 

max |/,M|= max II |sinc ^ | 



(fc) 



where /3 < 1 is a positive constant. And hence, when |f^| > tt/o^ 
(21) |/,(a;) - exp (-^) | < |A(u;)| +exp (-^) S.c^ p 
which implies ([20]) . 
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We next consider the case where |a;| < T^/a"^\ Taylor expansion shows 
that, when |a;| < ■K/a"^\ 



(22) In A(a;) = ^In j^sinc j j = " (p^ll' " ^ + S{u^] 

where 

II", l|4 , ,4 11^, ||6 , ,6 

^ ' 2880 181440 
is a uniformly convergent series on \uj I < Tr/af^- By ([22]), we now obtain 
that, when < T\'/a"^\ 

\2\ / ^ , ,2- 



,2 



Jk{i-}) = J__[sinc j = exp 1^ — j ■ exp |^ — — j ■ exp(-5(w)). 



Hence, 

(23) |/fc(u;)-exp 



24 



< exp(-5(a;))exp • l(exp ("^^ 

+ exp • |exp(-5(a;)) - 1|. 

Once, we prove that 

(24) exp • I exp (^-^^) " ^1 - l^^l ' M^k\\ 
and 

(25) exp(-^).|exp(-5(..))-l|<(i^^. 

Then, combining ([23]), ([M]) and ([25]), we obtain (|20]) . 
It remains to prove (I24p and (I25p . We first prove that 

V — 2F J ' '^^^ V — 2F ] ~ I ~ '^'=1 ■ l^'^l^'^ll- 

By Taylor's expansion, when lj^ < 24 • |ln|efc| j/cr^, 



when > 24 • |ln|efc||/cT^ 



This gives ([24 
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We next prove that (|25p . Note that, when < \/24:lnk/a, 
(26) exp (-^) |1 - exp(-5(a;))| < |1 - exp(-S(u;))| Sx,^ 
When \/241n/c/cr < |w| < 7r/ar\ we have 



a^uj^\ 1 



which imphes that 

(27) exp --- |exp(-5(a;)) - 1| < - < ^ ' 



24 y ' ' ' - ~ k 

Combining ([26]) and (f27l) . one derives (f25]) . □ 

Proof of Theorem (4]. Note that 



— y /fc(a;)exp(za;a;)da; = S(x|Hfc). 



Then 



°iax|y^exp -B(x|Hfc)| <^^|exp ) - fk{^)\duj 

.2, ,2^ 



/ ^ |exp T-^-^^ -/fe(w)|da; 



+ / ^ exp — - fk{u))\duj 

a, a, ' (It ' 

fc 1 1 

< ^ + k.l • llnk.ll • Infe + ^ ^ + le.l • |ln|6,|| • In A., 

where /3 = max{-,sinc ( f^ffey 1} < 1- Here, we use (]20p to obtain that 

/ ^ |/fc(c^) - exp (-^] \duj < + |e,| • |ln|e,|| • InL 

■/|a;|<^'"fe V 24 y ' A: 



Note that af'' /^/k, k = l,2,... , is a bounded sequence and 



, a^uj^, ^1 ^ ^ \/241nA; 

exp(-^)<^, fora.>^^. 
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Using a similar method as the proof of (f27|) , we have that 



/_, , - exp T-^-^^ldw < 

a —I I— (fc) 



fAn k ' 



To estimate 



/ |/fc(a;)-exp(-wV24)|f^a; < / \h{io)\duj+ [ exp(-wV24)dw, 
we use the facts of 



[ \fk{oj)\du;< [ 

J\uj\>-^ J\uj 



and 

. 1 



/ exp ( - 

Jh>^ V 
^1 



□ 



Theorem m impHes that the normahzed box sphne B^-l^m/) converges 
uniformly to a Gaussian function: 

Corollary 6. Suppose that 

Em,i = \l^_^,l/2,...,l/2]. 

m-i 2£ 



Then, for each fixed £, y^i — | • B{ym — converges uniformly to 

■ exp(— 6a;^), as m — )• oo. 
Proof. By the definition of box splines, we have 



Note that, for each fixed £, 11 ^"^'^ lln = 1. Then, Theorem H] shows that the 



box spline B{x\--^^=j^), and hence y"^ — | B{ym — ^x\E„^^£), converges 



uniformly to the Gaussian function y ^ exp(— 6x^). □ 

Remark 3. A well-known result is that \/rnBm{\/rnx) converges uniformly 
to y'^exp(— 6x^) with m — oo (see [HI [23]). In fact, the result can be 
considered as a particular case of Corollary [6| Note that 

y/m ■ B{y/mx\r.rnfl) = y/mBm{y/riix). 
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If we take ^ = in Corollary [Gj then we have that -y/m • B{y/mx\'Bmfi), and 
hence y/rnBm{y/mx) , converges uniformly to y'^exp(— 6x^) with m — )• oo. 

3.2. The asymptotic property of B-spline framelets. By changing 
variables, we can observe from Corollary [6] that B{x\'Em/) is close to 

6" 1 / 12-2;2 

: exp 



Recall that Theorem [3] says that 



Motivated by the two observations, we consider the relation between 'ip^^\x) 
and the ^th derivative of a Gaussian function G{x), which is defined as 

Gm/i^) ■= '^rn,e ' exp ' 



2m -£ 



where 



Let 



_ /6 



VT - i/2 ■ 4^ ■ 



(28) ) (x) := ^,GmA^ - ^-f), £ = l,...,m, 

where jm is given in ([3|) , and 

Then we have 

Theorem 7. Let m € N 6e given and 1 < i < m, and the framelet tp^/^^ 
be defined by its Fourier transform in ^ derived from B-spline of order m. 
Then, 

I ,(m). ^ ^{m), ^ (ln?n)^/2 
max max -i/;; ' [x) - G) [x)\ < -j- — . 

l<£<m xm ^ V ^1 ^3/2 

In order to prove Theorem [71 we need the following two lemmas: 

Lemma 8. The following inequality holds for every u, such that \u}\ > 
20 /i^. 



max 4 / „ \ ■ —7 ■ ■ exp — m — - 5- . . 



, . 1 
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Proof. For the convenient, denote 



and 



24 ■£ 



2m -£ 

For each fixed i € [l,m] n Z, the function is increasing on the inter- 
val [0, ti;^], while Fi is decreasing on [u£,oo). And hence Fi arrives at the 
maximum value at ui^. According to the inequality 

m\ / m • e N ^ 



we have 



Then, when 



, ^/ N ^ , 2(2m-£) 
2 6m 



21nm „ 2 

< £< -m, 



we have that 



In 16 - In 15 - "5 

1 ^ / N ^ , 2(2m -e) In f ^ 

In Feioji) < -In ^ < f^r Inm < -41nm. 

2 3m In 4| 

15 

This implies that, whenever 



21nm „ 2 

< i< -m 



In 16 - In 15 - "5 
holds, one has 

(29) maxF^w) < Fe{ue) < 

We next consider the case where ^ < £ < Using the inequality 



< 2"^, 



one gets that 



, , ln2 2e(2m-£) 
In FAujg) < m- - In— — -. 



Therefore, when ^ < £ < one has that 

1 N / ln2 £^ 2e{2m-£) , , 

inl'piLOe) < m In ^ —m ^ — 4inm, 

— 2 2 3-^ ^^"^ 

which implies that, whenever ^ < £ < the following holds: 

(30) maxF^w) < F^(w^) < 

ui m* 
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We now turn to the case |m < £ < m. For this case, we apply the inequahty 

\ / \ m—£ 



< 



m 



to obtain that 

, , m — £ , m • e £ , (4m — 2£) • e ^ 

In Fp(ujp) < In hi ^ < -m < -41n?n. 

- 2 m-£ 2 3e 

Hence, when < £ < m, we have that 

(31) maxFe{uj)<Fe{ooi)<^. 

We finally consider the case where 1 < £ < in i6- in 15 • Note that, when m is 
large enough, we have Ui < 20^J^^. When \uj\ > 20^^^, Fi{u)) reach the 

maximum value at cjq := "^^xf^^- Then a simple calculation shows that, 
when 1 < £ < 



2 Inm 



Inl6-lnl5' 

(32) max F^w) < F^(wo) < ^• 

U|>20,/i^ 



Combing 1^, jSOj), ((311) and ([32]), we conclude the proof. □ 
Lemma 9. The following inequality holds for every oj £M, 

V(T) , ,f , /LO\m-i ( (m -£/2)a;2\ , In^m 
— — ■ (j; • sinr, — sine — — exn — ■ — < 



max — — -„ — a; • smc — smc — — exp , 

\<l<m 4^ ''I \2) V4/ 24 

Proof. For the convenience, we only provide the proof for the case where 
a; > 0. The proof of the other case is similar. By Taylor expansion, when 
< a; < ^, we have 

smc(-j sinc(-j 
= 24 j^^H 2880 



Then, for 20 J <uj <^,we have 



uj\m~e _ /^x2^ / (m-£/2)uj' 



4. ^'-l^^^^y^) smc^-j -exp 2^ 



^ , f ( (m-£l2)u?\ / / (m - 7^/8)^^ 

^ ^^-rr- ■ ^ • exp ^ 1 - exp ' 



4^ V 24 / V V 2880 



i ( {m-£/2)u: 



,2X 1 



< 2^..^.exp -^-^^ <^,, 
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where the last inequahty is obtained by Lemma El Next, when < u) < 
20 ^i^, note that 

/ (m-7i/8)oj^\ / In^m 
1 - exp ( ) < 

and 



2880 / ~ m 



^, , , u) e ( {m-l/2)oj^ 



4^ "^l 24 



is a bounded function. Hence, for < w < 20^/^^, we have 



, uj ■ sine — smc — — exp , 
4^ I V2/ V4/ ^ \ 24 



Finally, We consider the case when > We assert that, when w > 
the following inequality holds: 



(33) max J , . |sinc - T"^ • |sinc - ^ < 



1/8 X"' 



i<£<myV^/ ^2/' ' V4/'~\37r 

With this assertion, we have 



-^■^ -Ismcj^-j sinc(^-J -exp' 



24 



< ^-0. -^Ismcy smc(-) | + | exp ^ j|J 

^ A/ « - smc - + -J • w -exp 



V 4 y I 4« I ^ V 24 



< J_ 



Here, the last inequality is followed by ([33]) and Lemma [8l It remains to 
prove ([33]) . Note that, when w > 

Isincf-'lr-^ • Isinc f-")!^ < ^ ^ < ^ ^ 

Then we only need prove 



max 



m\ 1 1 < / 8 • e^/^ 



i<e<m\ \ £ J (37r/4)'"-^ (37r/8)^ I 37r 
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/ \ m—i 

Applying the inequality < ( j , we have that, 



m 



1 1 /8-eV8^ 



max 4 / 1 • — — — J ■ — — — -T < 



i<i<m\j \i J (37r/4)'"-^ (37r/8)^ " I Svr / ' 
This proves that (|33p . □ 
Proof of Theorem [7]. Let 



Miuj) := max — — \uj\ - smc — smc — — exp —Im ) — - 

^ ^ i<e<m 4^ ' ' ' V2/ V4/ ^ 2^247 



and 



/ In TTh I In TTi 37r 

h := {u} e M. : U < 20\ }, /2 := G M : 20W < Iwl < — }, 

V m \ m 2 

Svr 

/3 := {w G R : > — }. 
Applying Lemma [9] and Lemma [HI we conclude that 



, , ^ (lnm)5/2 f 1 



respectively. By an argument similar to that leading to ()33p . we can obtain 
that there exists < 7 < 1 such that 



'h 

This leads to that 



/oo 
-00 

(Inm)^/^ 



M{uj)du}+ / M(tj)(itJ+ / M{uj)dLO 



< 



/2 -'/s "^^''^ 



□ 



Remark 4. It was proved in |17j that, for each fixed i, up to a normal- 
ization, a proper scaled V'^'"^ uniformly converges to the ^-order derivative 
of a scaled Gaussian function with m tending to infinity. Our result is in a 
different direction. In fact, we show that for sufficiently large m framelets 

ipi^^ ) • • • 5 i^i^'^ ) • • • ) V'm"^ uniformly in x and £ close to derivatives of consec- 
utive orders 1, . . . , m of a scaled Gaussian function whose scale depends on 
m. 
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4. Gaussian frame 

Theorem [7] leads us to consider whether a wavelet system generated by a 
finite number of consecutive derivatives of a properly scaled Gaussian func- 
tion forms a frame of -L2(M). In this section, we show that the frame property 
of X{'ij)^^\ . . . , V'm"^) can be transferred to that of X{G^^\ . . . , G^m^) where 

G^™^ is defined as follows: For each fixed m G N, we consider the following 
rescaled Gaussian function 

Gm/{x) = Cm,£ ■ exp 



12 • 

2m- i 



where 



and 



Crt 



vr y^m - i/2 ■ ¥ ' 



where jm is given in ([3|) , and 

G^'") = {gS'"\...,g(^)}. 

Before stating the following main theorem of this section, we recall the 
definitions of the frame and Bessel sequence. A family {fj}j£j C .^2(1^) is 
called a frame with bounds A and B if 

A\\ff<J2\{fJj)\'<B\\ff 

holds for all / S L2(K). li A = B, then {/j jjgj is called a A-tight frame. 
Moreover, a family C L2(M) is called a Bessel sequence with a 

bound R if 

Ek/'/.-)i'<^ii/ii' 

holds for all / e L2(M). 

Theorem 10. Let X{G^'^'^) be the wavelet system generated by functions 
ThenXiG^""^) is a frame system with frame bounds Am and Bm for 
sufficiently large m. Furthermore, the frame is close to be tight as m is 
sufficiently large. In fact, asymptotically, we have 

lim Am = lim Bm = 1. 

To prove Theorem [TOt we need the following theorem, which is proven 
in [16], together with several lemmas. 

Theorem 11 ( [E]). Let {fj}j^j be a frame of L2(^) with bounds A and B. 
Assume that {gj}j^j C -L2(M) is such that {fj — gj}j£j is a Bessel sequence 
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with a bound R < A. Then {gj}j^j is a frame with bounds ^ ^1 — 



and B 1^1+ § 
Let 

(34) <Af^=^f^-Gf\^ = l,...,m, and <I>M := M'"), . . . , 0^}. 

bmce is a tight frame with frame bound 1, to prove that 

is a frame, according to Theorem II l) we only need to show that X(<^(™)) is a 
Bessel sequence with a bound Rm — >• 0. An estimate of the Bessel bound of 
a given a sequence is provided in [21] that enables us to estimate the Bessel 
bound of (see also [H]). Let 

m 

(35) R^:= sup j;^^|0f)(2"o.)|-|,^f)(2"a; + 2A:^)|. 

i<|w|<2 i^fzznez e=i 

Then, for arbitrary / G L^(M), the following inequality holds 

(/,gX{<I>(")) 



i.e. Rm is the Bessel up bound of the system X(^>(™)). Next, we estimate 
Rm- For this, we need the following lemmas. 



Lemma 12. Let (^^™^ be the Fourier transform of (j)^^^ defined in ()34p . 
;n ti 

(i) 



T/ien i/ie following three estimates for \(p^^^\ holds: 



HI < \l{ J-^J^' M>20. 



(ii) 



(iii) 



max max \q)\ >, ■ 



l<£<m wgR " m 

Proof. First a simple calculation leads to 

\4>^r\^)\ = ivif^(-)-Gf)(^)i 
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For (i), when \uj\ > 20, a simple argument shows that 



exp -(m--)— < 



2' 24 



:^)^sinc(|r-^ 



sinc(:^)2^| < 



UJ 



which imphes that 



m 



2-m+^+l 



\U}\ 



For (ii), the Taylor expansion shows that, when |a;| < vr. 



In ( sinc( — ^sinc 



2' M' 
Then, when \uj\ < \/ — , 



21 



IsmcC-)"* 'sinc(- 



^ 2^24 ^ 8 ^2880 ^ ' 



exp(-(m-^)^) 



exp 



[m 1 



8 '2880 



<r 4 



which implies that 



Finally, the conclusion of (iii) can be obtained by Lemma [9] directly. 



m\ /ijj 



4 

m ■ Lo . 



Lemma 13. Let R„i be given by (I35p . Then 

In^ m 



Rm. S. 



and lim R„ 



m 



Proof. Let 
-Rm,i : 



■m,2 



m 

^<M<'2 kez\{o} n& e=i 



(2"a; + 2kTr) 



□ 



Then, we have that 

(36) Rra < Rm,l + Rm,2- 

To estimate Rm, we consider i?m,i and Rm,2, respectively. We first esti- 
mate i?m,,i- For this, we rewrite 

m 

Rm.,i = sup 5^J^|<^f)(2"a;)|2= sup + ^2(0;) + S3M], 



1< w <2 
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where 

m m. 

n>5 e=l - [log2 mj <n<5 ^=1 

m 

n<— [log2mJ £=1 

By (i) in Lemma \T2\ we obtain that, for 1 < |(^| < 2, 



= EEi'^f^(2"'-)i' 



n>5 £=1 



n>5 f=l ^ ^ 



< 



n>5 ^ ^ e=l ^ 

(- " 

Using (ii) in Lemma \T2\ when 1 < |f^| < 2, 

m m 

= E Ei^i"'(2"-)i^= E Ei^^^(^) 

n<— [log2'mJ (.=1 '^>U°g2™J 



< 



< 



^ ^ ( ^ ) (4 • 2") ■ ■ ( 24n 

n>[log2mJ ^=1 ^ ^ ^ 

4 \ 2 m / V , 



24n. / \^ £ y V4 • 2" 

n>[log2mJ ^ ^=1 

(4 \ 2 / 2\ ™ 

^ E (i) fi + (^ ' 

n>Llog2'TiJ ^ ^ \ ^ 

1 



< 



(2\ ™ 
l+(27^)j }n>Llog2mJ,mGZ+ 

being a bounded sequence and 



( 94n ) 



^ 1 256 1 ^ 1 



24n/ 256'^ 255 256Li°g2 H ' 
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Moreover, by (iii) in Lemma [T2| we have 

In^ m 



^3H = E E i<^^(2"^)i'^ 

£=1 -[log2mJ<n<5 

Combining the resuhs above, we obtain that 

In^ IT) 

(37) Rm,i= sup [S1M + 52M + 53H] < . 

l<|a;l<2 "T- 

We next turn to 



Rm,2= sup E EE(l'^^^(2-^)|-|'^^^(2"^ + 2^^)l)- 
To state conveniently, we set 



/3(2A;7r):= sup E E l<^l"'^(2"a;)| • |(A5'"^(2"w + 2A;7r)|. 

I<l^l<2^=l„g2 

Then 

Rm,2< E '^^2^^)- 

A;eZ\{0} 

Set ko := 10. When 1 < |A;| < A:o — 1, using the argument similar to the one 
in the estimation of Rm,i = /3(0), we can show that 

p(2kTT) <— for 1 < |A;| < feo - 1. 
m 

We claim that, when \k\ > ko, 

nm 

(38) /3(2A:7r) " 



< 



|2fc7rr/2- 
And hence, 

(\ 1 Qrrt 1 

l<|fe|<fco-l \k\>ko/ |fc|>fco I I 

Combing (f37|) and (f39l) . we obtain that 

m 

which implies the conclusion. 

Finally, we prove ([38]) . A simple observation is that /3(2A:7r) = f3{—2k7r). 
Hence, we only need consider the case where k > k^. For the convenience. 
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let 

m 

P+{2k7r) := sup ^ ^ (|<Af^(2-.;)| • |<^('")(2"u. + 2fc7r)|) , 

m 

P.{2k7r) := sup V V )(2"u;)| • |<^f )(2"a; + 2A;7r)|) . 

Then /3(2A;7r) < /3+(2A;7r) + /3_(2fc7r). To estimate /3+(2/c7r), we furthermore 
set 



/?+(2A;^) := sup V V ^2"^)! • |<A^^(2-^ + 2A;^)|, 

l<a;<2 



1 neZ+ 
m 



P+{2kTi) := sup ^ ^ (2"w)| •l.^r (2"'^ + 2A:^)|. 



— — 1=1 ng/+ 



Then /3+(2A:7r) = max{^|(2/c7r), /3+ (2A;7r)}. Note that k > ko = 10, by 
Lemma I12| 



/3| (2fc^) = sup ^ ^ 1 4>r (2"a;) I • 1 0f (2" w + 2A:7r) I 



i J {2'' + 2k7ry 



(2A;7r)'"/2 ' 
We next consider 



^ ) (2"^ + 2fc7r)'" 

. cm am 

< V ^ <V ^ - 

- _^ 2A;7r)™ - ^ 2™ • 2W2 . (2A;7r)™/2 



/3+(2A:7r) = sup ^ ^ |0P(2"a;)| • |(^P(2'^w + 2/fc7r) 



sup E ICH2"^)| • |<Ar(2"^ - 2A;7r)|. 

— — £=1 neZ-i- 



A simple observation is that max{|2"a;|, \ 2'^uj — 2k-K\} > kn > koir. Then 
using (i) in Lemma \T2\ we obtain that 

ir^(2"a;)|.|r)(2".-2^vr)|<(7).|^. 
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Set no := [log2(2fc7r)] + 5. Then, 



< 



sup J2 E l'^^^(2"c^)|-|<^f^(2"^-2A;vr)| 

l<<^<2 ^^-^ o<n<no 
<n<no ^=1 ^ ^ ^ ^ 0<n<no ^=1 



0<n<no e=l ^ ^ ^ ' 0<n<7io 
6™ 6"^ 



0<n.<no 



We next consider 



m 



sup E E I<^r(2"a;)|-|<^r(2"c^-2fe7r)|. 

<a;<2 (, , . 1 ^ 
— — £=1 no+l<n 



By (i) in Lemma [T2| when n > no + 1, 



l<a;<2 

Hence 



m 



sup 5] Yl I'^f^(2"^)|-|<^f^(2"c^-2fe7r)l 

- ^ ^ V £ y (2/c7r)'^ ' (2" - 2A;7r)™ 

- (2/c7r)'^ (2" - 2/c7r)™ ~ j ' " °' 



Therefore, 



This concludes that 



/?+(2A:7r) = max{/3+(2A;7r), /3-(2fc7r)} < 



(2A;7r)W2 " 

Using (ii) in Lemma [12] and a similar analysis with above, we can obtain 
that 

m 

/3_(2A;7r) = sup V V |<^{'^)(2"a;)| • |(A5'"^(2"w + 2A;7r)| 

l<|a^l<2^=l „<o 

< Tn • 

(2A;7r)'" 
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Table 1 . The numerical results of frame bounds of 



m 


2 


3 


4 


5 


6 


7 


8 


A 


0.3855 


0.5266 


0.5898 


0.6407 


0.6803 


0.7095 


0.7274 


B 


1.9020 


1.6239 


1.5179 


1.4390 


1.3811 


1.3403 


1.3159 



Putting everything together, we have that 

am Am om 

mkir) = p+ (2k7r) + /3_ (2fc7r) < ^ + m • , , , < . 

This proves (j38|) . 

□ 

Proof of Theorem IIOI Recall that 

and that is a tight frame with frame bound 1, where ^^"^i = 

Lemma US] shows that is a Bessel sequence with 

a bound Rm — )• 0. Then Theorem [11] leads that X(G*^™^) is a frame with 
frame bound Am = (1 — yJRmf' and = (1 + \/Rm)'^ as m sufficiently 
large. Furthermore, it can be close to a tight frame, since limm_j.oo A^. = 
limm^-oo Bm = 1 , which completes the proof. □ 

Remark 5. Although Theorem [10] confirms the case where m is sufficiently 
large, the result of Theorem 1101 seems to hold for small m [m can be as small 
as 2). For small m, combining (I35p and Theorem II H we can estimate the 
frame bounds of numerically. We list the frame bound estimation 

of 2 < m < 8, in Table 1, which clearly shows the frame property 

of for small m. For example, for m = 2, X(G*-^^) is a frame with 

frame bounds A ~ 0.3855 and B ~ 1.9020. Figure 2 shows that the graphs 

C9'J (2) 

of the functions of G{ and G2 , respectively. 

0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 

-0.6 




Figure 2. The graphs of GY' (x) = -y^32/7Tx expi-Ax"^) 
(left) and Gf\x) = ^27/(87r)(12x2 - l)exp(-6x2) (right). 
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